Magnetization transport and quantized spin conductance 
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We analyze transport of magnetization in insulating systems described by a spin Hamiltonian. 
The magnetization current through a quasi-one-dimensional magnetic wire of finite length suspended 
between two bulk magnets is determined by the spin conductance which remains finite in the ballistic 
limit due to contact resistance. For ferromagnetic systems, magnetization transport can be viewed 
as transmission of magnons, and the spin conductance depends on the temperature T. For antifer- 
romagnetic isotropic spin-1/2 chains, the spin conductance is quantized in units of order {giiB)^/h 
at T = 0. Magnetization currents produce an electric field and, hence, can be measured directly. 
For magnetization transport in electric fields, phenomena analogous to the Hall effect emerge. 

PACS numbers: 75.10. Jm,75.40.Gb 
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Transport of magnetization in various magnetic sys- 
tems has received considerable attention both theoreti- 
cally and experimentally Q, |^ IM Q . A spatially vary- 
ing magnetic field gives rise to a current of magnetic 
dipoles [3, , similar to the transport of electric charge 
driven by an electric gradient. Here we consider insu- 
lating magnets described by a spin Hamiltonian, where 
magnetization can be transported by excitations such as 
magnons and spinous without transport of charge. The- 
oretical work on such systems has been focused on the 
long-wavelength limit for magnets with translational in- 
variance 

In contrast, we propose to investigate magnetization 
transport in systems with broken translational invari- 
ance. In particular, we consider a quasi-one-dimensional 
system of finite length, e.g., a spin chain sandwiched be- 
tween two bulk magnets which act as reservoirs for mag- 
netization, where the magnetic field gradient is nonzero 
only over the system. Then, the magnetization current 
is determined by the spin conductance G which remains 
finite in the ballistic limit due to the contact resistance 
between the reservoirs and the system, in analogy to elec- 
tronic transport in mesoscopic systems 6]. This is in 
stark contrast to the spin conductivity which diverges in 
the ballistic limit due to translational invariance |3, 1^ ■ 
Here, we derive the spin conductance G for both ferro- 
magnetic (FM) and antiferromagnetic (AF) systems. We 
find that, for FM systems, magnetization transport can 
be viewed as transmission of magnons and the conduc- 
tance is temperature dependent. For the AF spin-1/2 
chain, the conductance has a value of order [g^.B)^ /h, 
where g is the gyromagnetic ratio and /i b the Bohr mag- 
neton. Further, spin currents produce an electric field 
which allows one to measure G. We discuss magnetiza- 
tion transport in an external electric field and show that 
phenomena analogous to the Hall effect exist. 

Ferromagnetic systems. - We first discuss a system 
with isotropic FM exchange interaction in a magnetic 
field B(xi) — BiGz- The spins occupy the sites of a 
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FIG. 1: (a) Proposed experimental setup for the measurement 
of a magnetization current /„. (b) A magnetic field difi'erence 
AS between the two bulk systems gives rise to 7m = GAB. 
(c) AB shifts the Bose functions ns(e) in the reservoirs Rl, 
R2. Magnons with energies e within the shaded region in Rl 
are not transmitted to R2. 



simple d-dimensional lattice with lattice constant a. 
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with J < 0. Here, is the spin operator of the spin with 
spin quantum number S at x.^, and (ij) denotes nearest 
neighbor sites. For spatially constant Bi = B > 0, the 
elementary excitations of the system are magnons with 
dispersion 
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Ek — g^J.BB + \ J\Sa k 
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which carry a magnetic moment —g^B^z- Here, k is the 
magnon wave vector. For temperature T <ti gfiBB/ks, 
the magnon density is small and the noninteracting- 
magnon theory is valid for all d. 

We now consider a setup for a magnetization trans- 
port measurement as sketched in Fig.^a). A spin chain 
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extends from x = — L/2 to L/2 and is suspended be- 
tween two large three-dimensional (3D) reservoirs, Rl 
and R2. L ^ a is sufBciently small that magnons propa- 
gate ballistically through the chain. The reservoirs nar- 
row adiabatically towards the chain ["transition region" 
in Fig. n^a)]. The system is still described by Eq. iQl, 
with the sites occupying a bounded region in space 
[Fig. n[a)]. A small spatially varying magnetic field 
5B{x)e^ with 5B{x) = -AB/2 (AB/2) for x < -L/2 
{x > L/2) is superimposed on the offset field Be^ for 
t > [Fig. mb)]. For < L/2, 5B{x) interpolates 
smoothly between the values ±A_B/2 in the reservoirs. 
The field gradient results in a magnetization current 
from Rl to R2. In linear response theory, can be 
expressed in terms of the spin conductivity a{x, x' ,10), 



Im{x,uj) = I dx' a{x,x' ,uj)dx' SB{x' ,uj). 
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To calculate Imix,Lu), knowledge of cr for x,x' £ 
[—L/2, L/2] is sufficient because dx'SB{x',uj) — in- 
side the reservoirs. For a quasi-one-dimensional system, 
due to the continuity equation, a is related to the sus- 
ceptibility X by (T{q,Lu) — —iujx{q,^) I (f' |1]- In the 
noninteracting-magnon approximation. 



{9^J■B) 



dk nsiek+q) - nBjek) 
2vr {ck+q - ek)/h + uj + iO' 



(4) 

Here, nsie) — l/[exp(/3e) — 1] is the Bose distribution 
function and P = l/ZcsT. In the limit w ^ of a 
dc field, from Eq. Q we find that linii^^o o'ix, x' ,uj) — 
(5M-b)^'^s(3Ms^)/^ is independent of x and x' . Inte- 
grating over x' in Eq. (P}, we find that 



nB{gfiBB)AB = GAB 
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is constant and depends only on the difference of mag- 
netic fields in the reservoirs, AB. Although magnetiza- 
tion is transported ballistically, the spin conductance G 
remains finite due to the contact resistance for magnons 
between reservoirs and the systern^similar to the related 
phenomenon in charge transport 6]. 

In FM systems, the magnetization current is carried by 
magnons. This allows us to reproduce Eq. (O from the 
Landauer-Biittiker approach y . The field difference AB 
switched on at t = results in a shift of the magnon en- 
ergies e in the reservoirs [Eq. (jSJ] and of the magnon dis- 
tribution functions (e) [Fig.^c)]. Hence, a nonequlib- 
rium situation is established. The magnetization in the 
reservoirs relaxes towards the new equilibrium values by 
magnetization transport from Rl to R2, i.e., the magne- 
tization current Im- AH magnons incident on the spin 
chain from R2 are transmitted into Rl In contrast, 
magnons with e e [gfiB{B - A_B/2), 5/15(5 -I- A_B/2)] 
are not transmitted from Rl to R2. This results in a net 



magnetization transport current 

Im = 9f^B / d€v{e)p{e)nB{€ + qhbB) 



-nB{g^iBB)AB - GAB, 



(6) 



where u(e) = dk^^k^/f^ is the magnon velocity and p{e) — 
\/hv{e) is the magnon density of states in the spin chain. 

If the system connecting Rl and R2 consists of sev- 
eral chains with finite interchain exchange J', G = 
(5M-b)^ Y.k^ nBigiJ-sB + ek^)/h, where Ck^ is the energy 
of the transverse magnon mode. At T = 0, G = 
because the system and the reservoirs are in the spin- 
polarized ground state. 

Antiferromagnetic systems. - As we show next, mag- 
netization transport in antiferromagnets is significantly 
different from fcrromagnets but similar to charge trans- 
port in Fermi liquids. In an AF chain with half- integer 
spin, the elementary excitations are massless, and we will 
show that G even at T = 0. The spin- 1/2 chain is 
believed to capture the essential features [icj, [ll|, [l^ ll3| . 
Thus, we now consider a spin-1/2 chain with isotropic 
AF exchange interaction J > in Eq. suspended be- 
tween two AF reservoirs For t > 0, a magnetic 
field B{x) is apphed along such that B{x) = —AB/2 
{AB/2) for X < -L/2 {x > L/2). By a Jordan- Wigner 
transformation and subsequent bosonization, the spin 
chain can be mapped onto a Luttinger liquid (spinlcss 
fcrmions). Then, at T = 0, the Euclidean Lagrangian 
flows into a massless free theory under renormalization 

group mm [3, 
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where K = 2^ v = {tt /2)Ja/h, and the homogeneous 
part of Sz is identified with c^s^/y^- The imaginary-time 
spin conductivity is (j{q, uJn) = {gfJ'B)'^{v /'KfiK)ujn/ (w^ -I- 
v^q^) 01 • However, in order to calculate G, it is not suffi- 
cient to evaluate the dc limit lo ^ Q oi (j{q, uj) because the 
elementary excitations change on propagation from the 
reservoirs (magnons) through the chain (spinous). Fol- 
lowing the related analysis for charge transport through a 
Luttinger liquid coupled to Fermi leads |3 j we model the 
transition from a 3D ordered AF state to the spin chain 
by spatially varying K(x) and v{x) in the Lagrangian 
Eq. (jZj). For simplicity, we assume that K{x) and v{x) 
change discontinuously from the values of the spin chain 
to the ones of a bulk antiferromagnet at a; = ±L/2 
[Fig.|2a)]. The values and Vb in the bulk region are 
choosen such that Eq. Q correctly reproduces the dy- 
namic susceptibility of a 3D AF ordered state. From the 
nonlinear sigma model description [l^ , we estimate v\j ~ 
V3Ja/n and Kb ~ 4\/3/7r. The spin conductance then 
follows from G = [(g^B)2/7r?i] lim^^^^o '^^^^^(a;, x', a;„) 



FIG. 2: (a) The transition from the 3D AF ordered bulk 
to the spin- 1/2 chain is modeled by spatially varying K{x) 
(solid line) and v{x) (dashed line) in Eq. Q. (b) setup for a 
transport measurement in an AF system. 



where the time ordered Green's function 

/"OO 

G^^{x,x\uj^)^ dTe-''^-^Trcl>{x,T)cf,{x',0)), (8) 

JQ 

must be evaluated for the inhomogeneous system in- 
cluding the transition regions |l6l |. For given x' € 
[— L/2,L/2], G^^{x,x' ,ujn) is obtained from the ansatz 
G^^{x,x' ,uJn) ~ acxp[ujnx/v{x)] + b eyi-p[—ujnX / v{x)] for 
the four regions x < — i/2, ~L/2 < x < x' , x' < x < 
L/2, and L/2 < x. The boundary conditions for the spin 
current are automatically satisfied by evaluating Eq. ||SJ). 
We find that limi^^^o iLinG^^{x, x' , ujn) — ^/2Ki, is inde- 
pendent of X, x' and of the parameters K and v of the 
spin chain. The spin conductance at T = 0, 



G 



hKh 



(9) 



depends only on the parameter Ki, of the bulk system. 

We next argue that Eq. @ remains valid also for a 
finite temperature and an offset magnetic field. Both 
analytical and recent numerical [3 analysis indicate 
that, even at finite T <C J/kB, Eq. Q describes the low 
energy behavior of the spin-1/2 chain correctly. Simi- 
larly, a finite offset magnetic field, g^sB <^ J, suppresses 
quantum fluctuations of Sz in the bulk, leading only to 
a slight decrease of Kb- Our result Eq. © is hence ro- 
bust both for finite temperatures and magnetic fields. 
In summary, an AF spin-1/2 chain suspended between 
AF reservoirs acts as a transport channel for magnetiza- 
tion with a spin conductance of order {gjiBY /h. For TV 
parallel spin chains with vanishing interchain exchange 
interaction, each chain acts as an independent transmis- 
sion channel and G increases by a factor N . The spin 
conductance is quantized in units of order {g^B)'^ /h. 

Preparation of a sample as shown in Fig. da) is chal- 
lenging. A promising strategy is the use of a bulk mate- 
rial with an intrachain exchange J much stronger than 
the interchain exchange, such as Sr2Cu03. If heated to 
temperatures T much larger than the Neel ordering tem- 
perature T/v, the spin chains decouple and magnetiza- 
tion is transported predominantly along the spin chains. 
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FIG. 3: (a) A current of magnetic dipole moment Im pro- 
duces an electric dipole field leading to a measurable voltage 
Vm- (b) Magnetic dipoles —g^LBez driven by a magnetic field 
gradient VB in an inhomogeneous electric field E(x) experi- 
ence a force F analogous to the Lorentz force. 



Hence, an AF wire heated to T > Tjv in its central part, 
but cooled to T ^ Tm at its ends [Fig. |2tb)] provides 
a realization of the system in Fig. ^a). Recent experi- 
ments [igj provide strong evidence that elementary ex- 
citations in various quasi-one-dimensional systems have 
mean-free paths of several hundred nanometers at tem- 
peratures up to 50 K. The mean- free path is limited by 
the defect concentration in the samples. For L < 1 /im, 
transport through the system shown in Fig. Efb) then is 
indeed ballistic as assumed above [Eq. 0]. 

Detection of spin currents. - A current of magnetic 
dipoles produces an electric dipole field. The electric field 
is most easily calculated by decomposing the magnetiza- 
tion current into contributions propagating at a certain 
velocity v, Im = gtJ>'B^T,n{v)v, where n{v) is the line 
density of magnetic dipoles with velocity v. For each 
the electric field in the laboratory frame is obtained by 
a Lorentz transform of the magnetic dipole field in the 
comoving frame. Summing over v, we find that the total 
electric dipole field [Fig. OJa)] 



Em(x) 



27r r2 



(0,cos2(/), — sin2(/)) 



(10) 



depends only on Im- Here, sine/) — y/r, coscj) — zjr^ 
and r = ^/y^ -V - For a numerical estimate, we now 
consider N parallel uncoupled AF spin-1/2 chains con- 
necting two AF reservoirs. With Eqs. © and (fTn|l . 

iE.„(x,i^<(»;f£^ = 4xio-.»^i. 

ZTT n. 4 r\mY rn 

(11) 

Even for moderate AB = 10 T and large r = 10 ^ m, 
the magnetization current transported by iV ~ 10** paral- 
lel spin chains leads to an electric field Em ^ 10~^ V/m. 
The voltage drop between the two points (0, r, 0) and 
(0,0, r) indicated in Fig. |2Ia) is then Vm = Em. r — 
10~i^ V, which is within experimental reach |20ll2l| . 
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Spin currents in electric fields. - A moving magnetic 
dipole moment also interacts with an external electric 
field E(x), leading to phenomena analogous to the Hall 
effect. A magnetic dipole —gfJ-B^z moving in an electric 
field acquires an Aharonov-Casher phase |22l | and the spin 
Hamiltonian is modified to 



+9l^B^BiSi^;,, 



(12) 



where s^ = Sj^x ± isj,y and 9ij ^ g^B J^' rfx • (E x 
ez)/hc^. Introducing magnon creation and annihilation 
operators, Eq. H12|l can be rewritten in terms of magnons 
with single-magnon Hamiltonian h. From Eq. (|12|) . 



(13) 



Here, we discuss only the classical motion of magnons 
propagating with velocity v = —Vx^x in a 2D system 
of finite width W in the y direction [Fig. |3Jb)], where 
Im — g^^B'nVxW , and n is the magnon density. From 
the equation of motion implied by Eq. (|13|) . one obtains 
the force acting on a magnon, F — —gfiB^[B — (v x E) • 
e^/c^] . The second term accounts for the interaction with 
the electric field. We now focus on E = E'{x,y,—2z) 
with E' = const. Then, the equation of motion of the 
magnons is formally identical to that of electrons in a 
constant magnetic field. Magnons are deflected into the 
By direction perpendicular to the transport direction e^;. 
Stationarity is reached when the magnon repulsion equals 
the driving force along Gj, due to the electric field. Taking 
into account only dipolar forces between the magnons, in 
the stationary state B — VxE'y/c^ is constant as function 
of y. The difference in magnetic fields AS = B{y = 
W/2) — B{y = —W/2) is related to the magnetization 
current density by the spin Hall conductance Gh, 



Irn „ AS 



gfiB'i^c^ A_B 



E' 



W 



(14) 



In the hydrodynamic regime, the drift velocity Vx is de- 
termined by the magnon scattering time r. At low tem- 
peratures, T is limited by impurities in the sample. For 
T on the order of 10^ - 10'^ ns, as measured for yttrium 
iron garnet (YIG) at 1-4 K (Ref.Il^), dxB = 10^ T/m, 
J = 200 K fcs, S" 1, and a = 1 A, the drift veloc- 
ity is Vx = 10 - 10^ m/s. A variation of electric field 
AE = E{y = W/2) - E{y = -1^/2) = 10^ V/m across 
the magnetic system then would lead to Ai? = 10^^ - 
10""'^ G resulting from the spin Hall effect. Thus, the spin 
Hall conductance Gh is within experimental reach. 
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